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A study of the adsorption of CO on late Ad and 5ci transition metal (111) surfaces (Ru, Rh, 
Pd, Ag, Os, Ir, and Pt) considering atop and hollow site adsorption is presented. The applied 
' functionals include the gradient corrected PBE and BLYP functional, and the corresponding hybrid 

' Hartree-Fock density functionals HSE and B3LYP. We find that PBE based hybrid functionals 

, (specifically HSE) yield, with the exception of Pt, the correct site order on all considered metals, 

04 ' but they also considerably overestimate the adsorption energies compared to experiment. On the 

other hand, the semi-local BLYP functional and the corresponding hybrid functional B3LYP yield 
very satisfactory adsorption energies and the correct adsorption site for all surfaces. We are thus 
faced with a Procrustean problem: the B3LYP and BLYP functionals seem to be the overall best 
choice for describing adsorption on metal surfaces, but they simultaneously fail to account well for 
the properties of the metal, vastly overestimating the equilibrium volume and underestimating the 
CO . atomization energies. 

Setting out from these observations, general conclusions are drawn on the relative merits and 
drawbacks of various semi-local and hybrid functionals. The discussion includes a revised version 
O . of the PBE functional specifically optimized for bulk properties and surface energies (PBEsol), a 

' revised version of the PBE functional specifically optimized to predict accurate adsorption energies 

(rPBE), as well as the aforementioned BLYP functional. We conclude that no semi-local functional 
is capable to describe all aspects properly, and including non-local exchange also only improves 
some, but worsens other properties. 

s 

'. ■ I. INTRODUCTION 

Q The adsorption of carbon monoxide has intrigued researchers in surface science for the last decades. This has obvious 

reasons. Carbon monoxide is a simple diatomic molecule, but despite its simplicity it shows a rich phase diagram on 
metal surfaces, [l], [3l The CO dissociation reaction involves a very simple reaction mechanism, but nevertheless the 
reaction barriers can be greatly varied by alloying or roughening the surface. It is fair to say that, together with 

^ ' hydrogen adsorption and dissociation, CO has become a classical benchmark for experimental and theoretical surface 
science techniques. This interest is also strongly driven by the importance of CO for many technologically relevant 
^O ' reactions.^ 

In this iight, it is unsatisfactory that state of the art density functional calculations fail to describe several aspects 
of the adsorption of CO on metal surfaces accurately. This concerns both the predicted adsorption site as well as 
the absolute magnitude of the adsorption energy. For the Cu, Rh and Pt (111) surface, the generalized gradient 
approximation (GGA) in the PBE version predicts that CO adsorbs at a high coordination site (typically the hollow 
site), whereas experiments unequivocally show that atop adsorption is preferred. For Ag and Au a degenerate site 
preference is predicted, which does not agree with experiments either. [J, Equally compelling is that adsorption 
energies are significantly overestimated using the PBE functional. [1, Q These discrepancies have first been noticed in 
a now classical paper by Feibelmann et al. and the name CO adsorption puzzle has been coined. [3| Since then there 
has been emerging evidence that the present local and semi-local functionals are not capable to describe the subtle 
d I balance between donation of charge to the substrate and back donation to the molecule correctly for many adsorption 
problems. However, it is still an open question whether this is a universal shortcoming of such functionals, or whether 
the functionals can be amended to improve the description of surface related properties without drastically worsening 
other important properties. [(| 

The common model to describe CO adsorption is the Blyholder model, which invokes interactions of the two 
CO frontier orbitals, the Scr HOMO (highest occupied molecular orbital) and the 2tt* LUMO (lowest unoccupied 
molecular orbital), with the metal states. Due to the interaction with the metal states bonding Scr- metal orbitals 
and anti-bonding 5cr-metal orbitals develop, and the latter are partly shifted above the Fermi-level of the metal, 
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TABLE I: Lattice constants a(A) of late Ad and 5d transition metals in the fee strueture using different functionals. Exper- 
imental lattice constants of Ru and Os have been calculated from the experimental hep volume assuming an fee structure. 
Experimental values are taken from Ref. [stI . Numbers in round brackets are the relative error (in percent) with respect to 
the experimental lattice constants. The experimental lattice constants have not been extrapolated to K, and the zero-point 
quantum fluctuations have not been included in the calculations. Inclusion of these contributions will change the values by 
approximately '^O.l % [36l ). 





EXP 


HSE 


PBE 


B3LYP 


BLYP 


Ru 


3.79 


3.76 (-0.8) 


3.80 (+0.3) 


3.82 (+0.8) 


3.87 (+2.1) 


Rh 


3.80 


3.78 (-0.5) 


3.82 (+0.5) 


3.85 (+1.3) 


3.90 (+2.6) 


Pd 


3.89 


3.93 (+1.0) 


3.94 i+l.S) 


3.99 (+2.6) 


4.04 (+3.8) 


Ag 


4.09 


4.14 (4-1.2) 


4.15 (+1.5) 


4.23 (+3.4) 


4.27 (+4.4) 


Os 


3.82 


3.82 (+0.0) 


3.85 (-f-0.8) 


3.87 (+1.3) 


3.92 (+2.6) 


Ir 


3.84 


3.84 (+0.0) 


3.88 (+1.0) 


3.90 (+1.5) 


3.95 (+2.9) 


Pt 


3.92 


3.93 (+0.2) 


3.97 (+1.3) 


4.02 (+2.5) 


4.06 (+3.6) 



causing a net bonding interaction (donation). Likewise, bonding 27r*-metal hybrid states become populated (back- 
donation). Simple symmetry reasons tell that the highly directional 5cr-metal interaction is particularl y st rong for 
atop adsorption, whereas the 2n* interaction is dominating for hollow site adsorption. [sl. [ol. Iiol. [Til. IT^. Il3l| 

Here we apply a set of now well established semi-local as well as hybrid functionals to the CO adsorption problem. 
Our focus is on the systematic variation across the periodic table, and on how hybrid functionals compare to the 
more traditional semi-local functionals. We show that none of the available functionals is capable to yield an equally 
satisfactory description of the metal and the CO adsorption. 



II. THEORETICAL AND COMPUTATIONAL METHOD 



The first principles density functional theory calculations and the hybrid functional calculations have been performed 
using the VASP code, PAW[14] potentials in the implementation of Kresse and Joubert[lB| and a cutoff energy of 
400 eV. The surfaces are modelled by a periodic four layer metal slab with c(2 x 4) symmetry, and the CO molecule 
is adsorbed on one side of the slab. We focus only on the atop and hollow fee adsorption sites, since the difference 
in calculated adsorption energies between the hep and fee hollow adsorption site is generally small, on the order of 
0.05 eV; see, e.g., Refs. [l^ [13, [iM E^- The Brillouin zone integration is performed on symmetry reduced F-centered 
6x6x1 grids (i.e. roughly 12 x 12 x 1 in the primitive cell). Since we want to concentrate on trends, all metals are 
considered in the fee structure, although Os and Ru crystallize in the hep lattice structure. 

The PBE[23 and BLYP [2ll [2211 functionals are used for the GGA calculations. For the hybrid functional calcu- 
lations, the HSE functional|23l. |24|. [25j has been applied. We use a variant of the HSE06 functional that observes 
the homogeneous electron gas limit and all important sum rules. Contrary to the conventional HSE06 functional, 
the screening parameter is set to w = 0.300 in both the semi-local as well as non-local part of the exchange 
functional, whereas the recommended choice is a; = 0.207 A^^. As has been shown previously, the specific choice of 
the screening parameter ui has very little influence on the total energies, but slightly affects band gaps.[2^ [2^, [l^l 
Throughout this paper we will use the acronym HSE for this functional. As we have shown elsewhere, for Cu, Rh and 
Pt the HSE functional yields the same site order and almost the same energetics as the more popular hybrid PBEh 
(sometimes also termed PBEO or PBEIPBE) functional Pl, ^3] (see Ref. HI for details). 

Furthermore the widely adopted B3LYP[21, 33] functional is applied in the present work. The difference to the pre- 
vious functional is predominantly in the correlation energy which is approximated using the semi-local Lee- Yang-Parr 
correlation. Systematic tests for selected extended systems are described elsewhere. [2^ [13, [H, [1^ H^l Implemen- 
tational details can be found in Ref. [2^. For the PBE and HSE functionals, the bulk metal lattice constants were 
optimized and used in the periodic slab calculations. In order to save compute time, the BLYP and B3LYP calculations 
were performed at the HSE geometries, which also agree best with experiment. 
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TABLE II: d band parameters for the 4d and 5d VIII B transition metals, is the center of the d band, is the d band filling, 
$ is the work function of the clean (111) surface, calculated using PBE. Numbers in round brackets refer to HSE calculations. 
$exp is the experimental work function. 41] "I>cxp of Ru and Os refers to the hcp(OOOl) surface. 





Ru 


Rh 


Pd 


td (eV) 


-2.59 (-3.17) 


-2.37 (-2.83) 


-1.84 (-2.09) 


rid 


6.50 


7.57 


8.71 


$(eV) 


5.1 (4.9) 


5.4 (5.0) 


5.4 (5.3) 


$cxp(eV) 


4.7 


5.0 


5.6 




Os 


Ir 


Pt 


€d (eV) 


-3.01 (-3.60) 


-2.93 (-3.39) 


-2.45 (-2.88) 


rid 


6.20 


7.24 


8.33 


$(eV) 


5.4 (5.4) 


5.5 (5.3) 


5.7 (5.6) 


$cxp(eV) 


4.8 


5.7 


5.7 



III. RESULTS 
A. Lattice constants 

The lattice constants of the considered metals are summarized in table HI We do not report the bulk moduli, since 
a one to one relation between them and the volumes exists (see e.g. Ref. [l^ [s^ for the PBE, HSE, and B3LYP 
case) : if the lattice constant is too small, the bulk modulus is overestimated, and vice versa. Hence the quality of the 
functional can be assessed by the predicted lattice constants. For Cu (Ref. ITsl. ITqI) . Rh, Ag, and Pt we can compare the 
PBE and BLYP lattice constants calculated using three different codes, e.g. VASP (this work), Quantum-ESPRESSO 
and Wien2k (Ref. H): in all cases, the agreement is very good typically within a few 0.1 %, except for Rh where the 
discrepancy between our PBE value and that in Ref. [y calculated using Quantum-ESPRESSO is 1.4 %. We relate 
this to the pseudopotential approximation applied in Ref. 

Table [J shows that the overall description is best on the HSE level, good on the PBE level, at best, modest using 
the B3LYP functional and quite bad using BLYP. The HSE functional always yields smaller lattice constants than 
the PBE functional, and the contraction is largest for half-filled d bands and smallest when the d band is entirely 
filled (Ag). Compared to experiment, this yields to a sizeable underestimation of the lattice constants for Ru, a 
sizeable overestimation for Ag and Pd, and fairly accurate values for the remaining elements. The B3LYP functional 
gives larger lattice constants with substantial errors for Ag, Pd and Pt (3.4 %, 2.6 % and 2.5 %). Finally, the BLYP 
functional yields the largest lattice constants with errors up to 4 % for Pd, Ag and Pt. Paier et al. have already shown 
in Ref. J35 that this is mainly related to significant errors in the LYP correlation functional for metallic systems. 

In summary, ad-mixing non-local exchange decreases the lattice constants, whereas the LYP correlation gives larger 
lattice constants than the PBE correlation. 



B. Metal d band parameters 

In table HIl we show the calculated metal d band parameters using PBE and HSE: the center of the d band, Cd and 
the metal work function, The experimental work function is also reported. The d band center can not be defined 
unambigously within a plane- wave approach. As we are mainly interested in trends {i.e. PBE vs. HSE), we choose 
to calculate the centre of gravity of the occupied d band (integration up to the Fermi level), which is a lower bound 
for the true d band center (see Ref. [l for details). In table [TTl we have also introduced a fractional d band filling (rid), 
calculated as the d charge contained in a sphere centered at the metal atom. The sphere radius has been fixed in 
order to contain exactly (within numerical error estimated to be 0.02 e) the number of valence s and d electrons in 
the atomic configuration. ^38<] Interestingly, it turns out that the volume of the sphere is almost equal to the volume 
per atom in the fee lattice (Wigner Seitz sphere), providing a physically sound basis for our definition of rid- 

We considered the bulk phase [iol. [i^ at equilibrium geometry using the PBE functional for calculating Ud. We 
will show below that this definition is helpful for discussing trends between and among the 4(i and 5d metals. From 
table ini we see that Cd calculated using HSE (e^^^) is always found at larger binding energies than Cd calculated 
using PBE (ef* ^^). This is related to the tendency of HSE to give a larger d bandwidth and a slightly weaker 
self-interaction within the d shell than PBE. (ISJ The common picture of reactivity of metal surfaces suggests that the 
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■-■ fee (PBE) 
•-• top (PBE) 




4d-metals 5d-metals 

3 I I I I I I I I 

Ru Rh Pd Ag Os Ir Pt 

atomie eonf. 4d'5s' 4d*5s' 4d'" 4d"'5s' 5d''6s" 5d'6s^ 5d'6s' 
6.6 7.6 8.7 9.6 6.2 7.2 8.3 

FIG. 1: Adsorption energies for the atop and fee hollow sites on late 4d and 5d metals for PBE and HSE. Experimental data 
with error bars are taken from Ref. 1441 The experimental adsorption site at low coverage is hollow for Pd and top for the other 
cases. 




4d-metals 5d-metals 

3 I I I I I I I I 

Ru Rh Pd Ag Os Ir Pt 

atomie eonf. 4d 5s 4d 5s 4d 4d 5s 5d 6s 5d 6s 5d 6s 
iij 6.6 7.6 8.7 9.6 6.2 7.2 8.3 

FIG. 2: Adsorption energies for the atop and fee hollow sites on late 4d and 5d metals for BLYP and B3LYP. Experimental 
data with error bars are taken from Ref. 44*. 



lower in energy Cd is, the smaller the adsorption energies should be. We will show below that this is in contradiction 
with our results for HSE. Although e^^"'^ < f^^^, the adsorption energies generally increase for HSE (see figure [1]). 

The HSE work function is generally smaller than the one calculated using PBE. This can be traced back to the 
induced surface dipole changes due to the inclusion of the Fock exchange. Indeed, hybrid functionals cause a complex 
redistribution of the charge density in the surface layer: an analysis of the difference between the PBE and HSE 
charge density (not shown here) shows that HSE functional depletes charge from the d^^ orbital in favor of d^z 
and dyz orbitals, thus decreasing the spill-out of electrons from the surface layer into the vacuum. This charge 
redistribution should also affect the initial steric repulsion, thus explaining why, for the top site, the HSE adsorption 
energies are larger than the PBE ones (see figures [1] [2] and [3]). The net effect of this complex charge redistribution 
is to reduce slightly the surface dipole and the work function. [43] The overall agreement with the experimental work 
function is good for both functionals. 



C. The gradient corrected PBE functional 

The central results of the present work are summarized in figures [T] and [21 where we show the adsorption energies 
for the atop and fee site, using PBE, HSE, BLYP and B3LYP. For B3LYP and BLYP, they are evaluated using the 
HSE equilibrium geometries. Experimental values are also shown with error bars. [441] We will first concentrate on the 
PBE results, which have been published in a similar form by Gajdos et al.^ For atop adsorption, the adsorption 
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FIG. 3: Adsorption energies for the atop and hollow site adsorption on late 4d and 5d metals shown versus the d band filling. 



energy progressively decreases towards the noble metals. The qualitative and to some extend quantitative behavior of 
the atop curve can be reproduced using the Hammer-N0rskov d band model. [39l. l45j in which the central parameter 
is the position of the d band with respect to the Fermi-level. It moves to lower energies as the d band becomes filled. 
A word of caution is in order here. Using our definition, the d band moves to higher energy as the d band filling 
increases (table |Tl|. This is not in contradiction to previous findings. ^39i] since the different behaviour is related to 
a different definition of the d band center. In Ref. [SOj the full d band is considered, while here we focus only on the 
occupied portion. The same trends as found in table [III have been reported in Ref. la. As we are mainly interested in 
trends, our following discussions are not affected by the particular definition. 

We note, however, that models based on the weak interaction limit, which forms the basis for the derivation of the 
Hammer-N0rskov approach, possess only limited value in the case of the strong interaction present for CO adsorbed 
on a transition metal surface, in particular, at the hollow site.Q] In fact, the d band model describes well the overall 
trends in the periodic table, but it does not account for the site dependent geometric effects (i.e. top vs. hollow). A 
convenient way to present the results for the 4c? and 5d transition metals is to show all values in a single graph and 
to choose as the abscissa the occupation of the d band (figure [3]) according to the previous definition. Interestingly 
for the top site, all adsorption energies fall on a single curve, suggesting that for top site adsorption the filling of the 
d band and not the position of the d band is the most sensible parameter (descriptor) to characterize the interaction 
strength. Figure [3] clearly shows a decreasing preference for atop binding in the sequence Os<Ru~Ir<Rh<Pt<Pd. 
These findings agree with experimental evidence within experimental errors (see Ref. 44 for further details) and other 
theoretical calculations. |46i] Furthermore, for many properties, for instance catalytic properties, Pt behaves similar to 
a fictitious element between Rh and Pd, as the d band filling is roughly in-between the latter two elements. According 
to conventional wisdom, the 5a interaction dominates for atop coordination. Hence, the decrease of the adsorption 
energy with increasing na filling can be explained in terms of occupation of antibonding 5a — d metal states that 
destabilize the molecular surface bond. 

Quantum chemical calculations have shown that the back donation mechanism is the most dominant for the hollow 
site (27r* — d interaction). [U 111] We find, in agreement with previous studies, 0] that the hollow site adsorption energy 
generally increases towards the late transition metals, with a final sharp decrease for the noble metals. This decrease 
for the noble metals is related to a fully occupied d band, located several eV below the Fermi-level. Therefore, for noble 
metals it becomes difficult to shift antibonding 5a — d states above the Fermi-level, and additionally the interaction 
between the deep and contracted d orbitals and the 27r* frontier orbitals becomes small. It turns out that sensible 
descriptors are the work function and the d band filling. In fact, previous studies have shown that the back-donation 
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interaction is very sensitive to the relative position of the 27r* orbitals and the d band, which also depends on the 
metal work function $.[4^, [1^ Clearly, a larger $ value at fixed or almost constant d filling {i.e. Os vs Ru; Ir vs 
Rh; Pt vs Pd) implies a downward shift of the Fermi energy, decreasing the back-donation into the 2n* — d bonding 
orbitals, thereby decreasing the fee adsorption energy. As a result, the CO adsorption at the fee site is less stable 
on Os, Ir and Pt than on Ru, Rh and Pd. When increasing the d filling, that is moving from left to right along 
the periodic table, the back-donation is also strengthened: since the bonding 2tt* — d orbitals become increasingly 
occupied, the molecular surface bond is stabilized. Counteracting this trend is the change in the work function, which 
increases towards the late transition metals. Finally, steric repulsion, which is always stronger at the hollow site than 
at the atop site, also plays a role. It decreases when moving to the right of the periodic row.[5Q| Although we found it 
impossible to disentangle the aforementioned factors, the calculations show that the combined effect of these trends 
is that hollow site adsorption becomes stronger towards the late transition metals. Last, but not least, we want to 
note that the diff'erent qualitative trends for the atop and fee adsorption energies as a function of the d filling reflect 
the different dominant interaction term at the atop and fee site, that is donation and back-donation, respectively. 

We now return to figures[T]and [21 It is clearly born out that the top site is preferred for early transition metals, but 
moving along the transition metal series the hollow site becomes progressively more stable. For the PBE functional, 
a change in site preference is observed between Ru and Rh for the Ad metals, and Ir and Pt for the 5d metals, earlier 
than observed in experiment. The noble metals behave somewhat differently. As already argued above, this is related 
to a very weak 5cr — *■ d donation and weak d — > 27r* back-donation. Here bonding is dominated by the s and p orbitals, 
but again the PBE functional fails to predict the right site order. A 0.1 eV shift of the fee and top curves against 
each other would give the correct site preference for all metals, e.g. top for Rh and Pt, but maintaining hollow site 
adsorption for Pd 

D. The hybrid HSE functional 

The hybrid HSE functional shifts the 27r* orbitals towards the vacuum level and the 5a orbital towards stronger 
binding energies. jll] In combination with a downshift of the d band center, caused by a reduced self-interaction 
within the d shell due to the inclusion of Fock-exchange, one might naively believe that the interaction between the 
2tt* orbitals and the metal d states should weaken for hybrid functionals. Although this conjecture is confirmed for 
the noble metals, for transition metals the weaker interaction is counteracted by the much increased d band width of 
the metal and a concomitant decrease of the work- function from PBE to the HSE functional (see table |TT| . 

It is clear from both, figures [U [5] and [31 that the adsorption energy increases for most adsorption sites using HSE: 
(i) the increase is pronounced and roughly constant for the top site on the transition metals; (ii) for the hollow site a 
slight destabilization is visible for Os and Ru, and a slight stabilization is visible for Pd and Pt. It is remarkable how 
smooth the changes along the transition metal series are, supporting the high precision of the present calculations 
and excluding numerical artifacts. As a result of these shifts, the crossing points between the hollow and top site 
move towards the right, for both Ad and 5d metals, but not sufficiently so to favor atop adsorption on Pt(lll). For 
the other critical cases, i.e. Rh, Ag (and Cu, not included here but discussed in Ref. |18i), HSE is capable to recover 
the correct site order, but only just. The main reason why HSE fails to predict CO adsorption above the top Pt site 
is that the hollow site also becomes more stable for the late transition metals, in particular Pd and Pt, whereas the 
increase in stability is small at the hollow site for earlier transition metals. This increased stability may be related 
to the tendency of HSE to increase surface energies (see below) and the tendency to localize the charge at the metal 
atoms, causing charge depletion at the hollow sites. In figure [H we show the difference between the PBE and HSE 
charge density, plotted for a (111) plane ^ 0.7 A above the metal atoms of the clean surface, that is, roughly in 
between the CO molecule and the metal surface layer. The figure indeed indicates charge depletion above the hollow 
site (red regions), suggesting a decrease of the steric repulsion using HSE. Furthermore, figure [H suggests that, going 
from PBE to HSE, the hep site should increase more in stability than the fee site. This is confirmed by our previous 
calculations for Rhls^ that gave a stabihzation of AEf^^^^'^^'^ ^ 0.007 eV for the fee site, and AE^^^^^'^^'^ = 0.027 
eV for the hep site. 

The most unsatisfactory result is that the adsorption energies increase substantially for the transition metals from 
PBE to HSE, which leads to fairly strong disagreement with experiment on the HSE level — even more so than for 
the PBE functional. We note that the same site preference, but even larger adsorption energies are obtained using 
the standard PBEh hybrid functional (also termed PBEO or PBEIPBE functional). This functional also includes 
25 % non-local exchange but, contrary to the HSE functional, does not truncate the long range part of the non-local 
exchange. [19] Since the only difference between the PBE and PBEh functional is that the latter one replaces 25 % 
of the local approximation to the exchange by the Hartree-Foek exchange, the conclusion is that inclusion of 25 % 
non-local exchange, bare or long-range screened, hardly improves the description of CO adsorption on metal surfaces. 

This confirms our previous study, [ll, [l^ but contradicts a recent work, where it is claimed that the CO adsorption 
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FIG. 4: Difference between the PBE and HSE charge density for Rh and Pt clean surfaces plotted in a (111) plane ~ 0.7 A above 
the surface plane (i.e. between the C atom of CO adsorbed at the hollow site and the surface metal atoms). Blue (red) color 
corresponds to regions where the PBE charge density is larger (smaller) than the HSE one. Three isolines are drawn to guide 
the eyes, "fee" and "hep" indicates the location of the fee and hep hollow site, respectively. 

puzzle has been solved for Pt using hybrid functionals. [Slj In order to explain the contradiction, we have calculated 
the adsorption energies using norm-conserving and PAW potentials and we present the results in the Appendix. We 
show that the conclusions of Ref. [Sll derive from an artifact of the pseudopotential approximation, and all-electron 
methods still yield the wrong site order for CO on Pt(lll) using the PBEh functional. 

E. The hybrid B3LYP functional 

Our results also somewhat disagree with other reports, which have applied the B3LYP functional and found stronger 
preference for atop adsorption and weaker adsorption energies on Cu and Pt.[l3, [H, [H, [s^l To resolve this issue, 
B3LYP values are reported in figure [H In fact, the B3LYP functional shifts the crossover point even further to 
the right and lowers the adsorption energies significantly compared to the HSE and PBE functional. Since these 
"improvements" can not be related to the inclusion of non-local exchange (see above) , the difference must stem from 
the semi- local part, and in all likelihood the Lee Yang Parr correlation energy (see e.g. Ref. Issh . In agreement with 
this conjecture, the semi-local BLYP functional also predicts the right site order on all considered metals. In summary, 
improvements from the PBE to the B3LYP functional are mostly related to the LYP correlation functional, and only 
partly to the inclusion of non-local exchange, confirming the recent work by Alaei et a/.[g] 

The seeming "improvements" for the BLYP and B3LYP functional over PBE and HSE, however, come at a consid- 
erable price best illustrated in table HI The equilibrium volume of Ag is overestimated by 12 %, that of Pt and Pd by 
about 9 %, an error that is by no means acceptable. Somewhat less troublesome but still worthwhile mentioning is 
that adsorption on the noble metals becomes slightly endothermic using B3LYP and BLYP (no geometry optimization 
was performed for the B3LYP case). We do not want to overemphasize this point, since presently DFT functionals 
generally lack weak van der Waals like dispersion forces, which might contribute to the CO bonding on noble metals. 
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Rh Pt 




FIG. 5: Surface energies a, adsorption energies Eads for atop and fee hollow sites, and energy difference between top and fee 
sites ABtop-fcc for CO on Rh and Pt(lll) versus theoretical lattice constant a (A) calculated using a set of different functionals. 
The experimental surface and adsorption energ y is taken from Ref. [s^ and Ref. lU respectively. The experimental AJ5top-fcc 
for CO on Pt(lll) has been inferred from Refs. Is^. IstI Surface calculations were performed at the theoretical lattice constant 
corresponding to the applied functional. 

IV. DISCUSSION 

It has recently been suggested that the quest for improved semi-local functionals is not yet over, and that it might 
be possible to design semi- local functionals that better describe the adsorption on surfaces. [6] Unfortunately this is a 
misconception, more precisely, for semi-local functionals, it seems impossible to improve the description of adsorption 
without worsening other properties. To discuss this issue we have extended our study to three additional functionals, 
the rPBE functional by Hammer and N0rskov,[58'| and the recently suggested PBEsol[6^ and AMOS'so'l functional. 
The last functional is based on the subsystem approach and designed to yield accurate surface energies. The other 
two functionals are modifications of the standard PBE functional. These modifications are in the opposite direction: 
the PBEsol functional follows similar ideas as the AM05 functional in order to improve surface energies, whereas the 
rPBE functional improves formation energies and adsorption energies. 

Most gradient corrected functionals essentially modify the LDA exchange and correlation energy to account for 
areas of large gradients. This is usually done by increasing the exchange and correlation energies locally in those areas 
where the gradient is large. For the "canonical" PBE functional, the enhancement factors have been determined using 
sum-rules, and the result is a well balanced functional that describes molecular and solid state properties equally well. 
There are however failures, in particular small molecules are over-bound, whereas solids involving elements beyond 
the second row are usually under-bound. Resultantly, for the metals considered here the lattice constants are too 
large and the surface energies much too small (underbinding) . To highlight this, figure [5] plots the surface energy for 
all considered functionals versus the theoretically predicted lattice constants of Rh and Pt. Similar results are to be 
expected for the other transition metals. 
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The AM05 functional has been designed to describe jeUium and jeUium surface energies accurately. The PBEsol 
functional mimics the AMOS functional by reducing the enhancement factor in the area of strongly varying gradients, 
thus approaching the LDA description. [60| As a rule of thumb, both increase the cost for surfaces, and hence systems 
involving surface areas become less favourable. Atoms are destabilized, the formation energies of molecules increase 
(molecules possess a smaller surface area per atom than atoms) and the atomization energies of bulk materials become 
much larger (no surface area). Furthermore bulk lattice constants shrink to reduce "internal" surfaces. The AM05 
and PBEsol functionals shift the focus towards a proper description of solids, but in doing so small molecules are 
strongly overbound compared to atoms using the PBEsol and AM05 functional. [60] In figure [5] (a), the AMOS and 
PBEsol functionals yield smaller lattice constants and larger surface energies than the PBE functional; both are now 
in quite reasonable agreement with experiment. 

The rPBE functional does exactly the opposite: it increases the enhancement factor in the area of strongly varying 
gradients, and it therefore moves further away from the LDA description. As a rule of thumb, the rPBE functional 
decreases the cost for surfaces. Therefore, molecules are destabilized compared to atoms, bulk materials become less 
stable than molecules, and bulk lattice constants increase. The decrease in the surface energy and increase in the 
lattice constant is again visible in figure [5] (a). The BLYP functional, characterized by an even larger enhancement 
factor (see e.g. Ref. |6|), continues this trend. 

We can now understand, why the adsorption energy decreases from AMOS and PBEsol, over PBE, rPBE to BLYP 
[Figure [5] (b)]: when CO adsorbs on the surface, the total surface of the combined system CO + surface is effectively 
smaller than that of the isolated CO and the bare surface.^] In AMOS and PBEsol, surfaces are expensive, and thus 
adsorption processes, which always lead to a reduction of the effective surface area, yield a larger energy gain. For 
rPBE and even more so for BLYP, surfaces are cheap, and less energy is gained during adsorption processes. 

The change in the energy difference between the top and the hollow site [figure[n](c)] can be understood by reasoning 
along the same line. Atop adsorption reduces the surface area very little essentially along a single bond between the 
CO and the surface, whereas hollow site adsorption reduces the surface area along three bonds between the three 
neighbouring surface atoms and the CO. If surfaces are "expensive" (PBEsol), hollow site adsorption is preferred, 
since hollow site adsorption reduces the effective surface area most. If surfaces are "cheap", top site adsorption is 
preferred. 

The binding energies, surface energies, lattice constants, and adsorption energies are thus all linearly related, and 
we could plot any quantity versus the other and we will always observe a linear relationship. In figure O we have 
chosen the lattice constant as independent variable, but the surface energy would work equally well (and arguably 
might be the better control parameter). Within the family of gradient corrected functionals, PBE seems to remain 
the best compromise. The AMOS and PBEsol functionals describe the lattice constants of the solids best and also 
yield improved surface energies, but unfortunately they overestimate the heats of formation, and important for the 
present study, the binding energies of molecules on surfaces. The BLYP functional arguably yields the best adsorption 
energies and the right site order, but for the price of crossly wrong bulk energies, much too large lattice constants, 
and much too small surface energies. 

The hybrid functionals constitute a deviation — albeit only a subtle one — from this general trend. The inclusion of 
exact exchange also clearly increases the cost for surfaces. Resultantly, HSE leads to a smaller lattice constant than 
PBE, and B3LYP to a smaller lattice constant than BLYP. Even the average adsorption energies follow roughly the 
same trend as for semi- local functionals: in particular the adsorption energy for the top site lies on the same curve as for 
the semi-local functionals. The only major deviation is visible for the energy difference between the top and hollow site. 
This we believe to be related to the shift of the 2tt* LUMO orbital towards the vacuum level, decreasing the capability 
of the 27r* orbital to accept charge from the substrate. As we have argued before this generally decreases the adsorption 
energy, with a stronger effect for the hollow site than the for top site. On passing, we note that on the basis of surface 
energy arguments alone, one would expect hybrid functionals to yield similar atomization energies for molecules and 
solids as a semi-local functional predicting similar lattice constants. This would imply that the atomization energies of 
molecules were significantly overestimated and that of solids slightly overestimated for the HSE functional. However, 
as we have reported previously, this is not the case, i.e. the atomization energies of molecules agree quite well with 
experiment, [34j whereas those of solids arc underestimated. |,26„ ,3S.] The reason for this is that the inclusion of the 
exact exchange stabilizes spin polarized solutions compared to non-spin polarized solutions, e.g. spin polarized atoms 
are much more stable for hybrid functionals than for semi- local exchange correlation functionals. For elements and 
d metals this leads to a slight and sizeable underestimation of the atomization energies, respectively. |26|] Neglecting 
these important spin-polarization effects, hybrid functionals do behave reasonably similar to semi-local functionals, 
with the tendency of increasing the cost for surfaces compared to the corresponding semi-local functionals. 

Clearly, however, none of the functionals can predict the surface energies and the adsorption energies well at the 
same time. From our point of view, the good prediction of the B3LYP and BLYP for adsorption energies and the 
adsorption site is not related to a better description of the underlying physics but rather accidental and fortuitous. 
By reducing the surface energy to an, in fact, unrealistically low value — worsening the lattice constants — the right 
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site order is eventually recovered, but there is little doubt, the reduction of the surface energy is not in accord with 
experiment. We believe that the main error of semi- local functionals is in the description of the strength of the back 
donation to the 2tt* orbital which tends to be overestimated using any of the semi-local density functionals. Hybrid 
functionals improve on that aspect, but unfortunately simultaneously increase the d band width and concomitantly 
the overall adsorption energies. 



Overall, the situation is unsatisfactory. Considering CO adsorption alone, the B3LYP and BLYP functionals seem 
to offer a very good description, but the price to pay is that the description of the metals, in particular d metals, 
is unsatisfactory. Errors of 10 % in the volume are unacceptable by today's standards, and worse, the atomization 
energies of d metals are wrong by up to 50 % for the B3LYP functional. [Ill Previously published seemingly good 
B3LYP results [l^ilHjIll, [131 must be reconsidered, when the broader picture is taken into account. B3LYP and BLYP 
fail for metals and surface energies: we feel there is hardly a point to predict surface properties, if the underlying 
substrate and surface is not accurately and properly described. 

The HSE functional performs significantly better for metals, yielding a better overall description of the equilibrium 
lattice constants than most semi-local functionals. But again, atomization energies for transition metals are wrong 
by typically 20 %, 26] and the overall description of CO on transition metal surfaces improves only little compared 
to the PBE functional. In particular, the increased adsorption energies are unsatisfactory. Overall, we are forced 
to conclude that for metals and metal surfaces, hybrid functionals are hardly a major step forward. In accord with 
our previous conclusions, [l^ we believe that the inclusion of a sizable amount of non-local exchange is not capturing 
the physics properly in metals. The non-local exchange interaction is strongly screened in metals, both in the long 
range as well as medium range. In combination with a semi-local correlation functional, the hybrid IISE06 functional 
overestimates the exchange interactions at medium distances, [isl [26| 

Finally, we have shown that trends from one to the other semi-local functional can be easily understood considering 
the cost for creating a surface. For the recently suggested AMOS and PBEsol functionals, surfaces are, in agreement 
with experiment, fairly expensive favouring bulk-like behaviour, whereas rPBE and BLYP underestimate surface 
energies and shift the stability towards smaller fragments (molecules and atoms). It is really puzzling that AM05 and 
PBEsol yield good bulk lattice constants and surface energies but overestimate adsorption energies, whereas rPBE 
and BLYP underestimate binding in solids and surface energies but describe adsorption energies well. A functional 
that could resolve both issues would be a tremendous step forward but we firmly believe that semi-local functionals 
can not achieve this goal. 



This work was supported by the Austrian Fonds zur Forderung der wissenschaftlichen Forschunq. A.S. thanks J. 
Paier for interesting and fruitful discussions. Figure [3] has been done using the XCrySDen package. [61| 



VII. APPENDIX. EFFECT OF PSEUDOPOTENTIAL APPROXIMATION IN HYBRID-FUNCTIONAL 

CALCULATIONS 

In a recent work[5l| it was suggested that the proper site preference for CO on Pt(lll) can be obtained using the 
PBEh hybrid functional. This observation contradicts our own finding. Il8| and we show here that the previous result 
is an artefact related to the pseudopotenti al ap proximation applied in Ref. ISll To understand the discrepancy between 
our own work and that published in Ref. [5l|, we generated a norm conserving pseudopotential (NC-PP) mimicking 
the FHI NC-PP applied in Ref. Isil. Figure [6] shows the density of states (DOS) evaluated using this normconserving 
pseudopotential and the PAW method. In our implementation, the PAW method has several advantages over NC-PPs. 
First, the exact shape of the all-electron orbitals is restored avoiding any shape approximations in the charge density, 
potentials and non-local exchange interaction. Second, although the PAW method applied here is using frozen core 
orbitals imported from a GGA calculation for the atom, the exchange and correlation energy is re-evaluated exactly 
and consistently using all electrons (core and valence), whereas in the pseudopotential approximation the hybrid 
functional is applied to the valence electrons only, and the core-valence interaction is implicitly calculated on the 
GGA level. Clearly, both approaches result in a different density of states for bulk Pt using the PBEh functional (see 
figure [HI). We note that no difference is observed on the GGA level, for which the potentials were initially generated. 
Furthermore, we observe that the PAW results are entirely robust and do not change upon treating the 5s or 5p shell 
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Energy (eV) 

FIG. 6: Electronic density of states (DOS) of Pt calculated using a normconserving (NC) potential, a normconserving potential 
treating the 5s and 5p as valence (NC— 5sp) and using the PAW method for the PBEh functional. The Brillouin zone was 
sampled using 8x8x8 k-points and a smearing width of a = 0.3 eV was applied. 



TABLE III; Adsorption energy of CO on Pt(lll) evaluated for the top and fee site using a normconserving (NC) potential, 
a normconserving potential treating the 5s and 5p as valence (NC— 5sp) and using the PAW method. The energy cutoff has 
been fixed to 1000, 1100 and 600 eV respectively. Numbers in round brackets are taken from Ref. 51 . 





NC 




NC-5sp 


PAW 


PBEh top 


-2.07 (- 


1.80) 


-1.92 


-1.88 


PBEh fee 


-1.85 (- 


1.67) 


-1.89 


-1.88 


PBE top 


-1.56 (- 


1.53) 


-1.56 


-1.54 


PBE fee 


-1.70 (- 


1.64) 


-1.67 


-1.66 



as valence (not shown), whereas placing the 5s and 5p shell in the valence for the NC-PP (NC— 5sp in figure[6]) yields 
results that are practically identical to the PAW results. The implicit treatment of the core- valence interaction on 
the PBE level for the standard NC-PP causes fairly large errors. 

In the second step we evaluated the CO adsorption energy at the fee and top site for a -\/3 x \/3 slab using these 
three potentials and 4x4 k-points (similar technical parameters were applied in Ref. Isil) . The results are summarized 
in table IIIII At the PBE level, all three potentials give very similar results to those of Ref. [Sll. Large differences 
are observed at the PBEh level. The standard NC potential yields the same site order as in Ref. [Sl] (but different 
adsorption energies), whereas the NC— 5sp potential — treating the 5s and 5p electrons explicitly as valence — and 
the PAW method show that the top and fee site are practically degenerate for this setup. For a lower coverage using a 
c(2 X 4) surface cell and using more k points, the fee site is consistently preferred by 0.04 eV for the PBEh functional. 
We conclude: the PBEh functional does not predict the right site order on Pt(lll) in the low coverage limit. Also 
the adsorption energy is significantly overestimated by the hybrid PBEh functional, and overall the description is 
thus not improved over the GGA case. [3 Furthermore, normconserving transition metal potentials generated using 
LDA/GGA functionals must be used with great care in hybrid functional calculations. We believe that this rule also 
applies to elements beyond the third row (4p, 5p and 6p elements), where the core valence overlap is often appreciable. 
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